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XXX. On Simultaneous Differential Equations of the First Order in which the Number of 
the Variables exceeds by more than one the Number of the Equations. By Geokge 
Boole, F.R.S., Professor of Mathematics in Queen's College, Cork. 

Received June 19, — Eead June 19, 1862. 

It is a fundamental proposition of analysis that a system of n differential equations of 
the first order containing %+l variables admits of n integrals, each of which is 
expressed by a function of the variables equated to an arbitrary constant. 

But when a system of n differential equations of the first order connects n-\-r vari- 
ables, r being greater than unity, no existing theory assigns in a general manner the 
number of theoretically possible integrals of the above species, or shows us how to dis- 
cover them. Yet such cases are of great importance. 

I wish to develope here the theory of a method for the solution of the above classes 
of equations, which was published by me in the 'Proceedings of the Eoyal Society' for 
March 6th of the present year, and which enables us to assign the number of theoreti- 
cally possible integrals, and to reduce their discovery to the solution of a system of 
simultaneous differential equations equal in number to the number of integrals, and 
expressible as exact differential equations. 

The solution of the problem as thus reduced may be effected by known methods, but 
I have thought it desirable to discuss this part of the subject also in direct sequence to 
the other, and in conformity with its method. 

Of the Connexion between ordinary and partial Differential Equations. 

It has been found convenient, in researches bearing upon the general theory of dif- 
ferential equations, to use the term 'integral' in two distinct senses, viz. to denote, as 
above, a relation satisfying the differential equation or system of equations, and expressed 
by the equating of a function of the variables to a constant, and to denote the function 
itself. The particular sense intended will always be shown by the connexion. 

With this convention two systems of differential equations will be said to be equi- 
valent when they have in either of the above senses (and the one implies the other) 
the same system of integrals. This will explain the meaning of the following pro- 
position. 

Pboposition I. — A system of n ordinary differential equations of the first order con- 
necting n+r variables may be converted into an equivalent system of r linear partial 
differential equations of the first order. 

Let x„ x 2 ,... x n + r be the variables, then the supposed given system of differential equa- 
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tions may, by algebraic solution with respect to dx„ dx 2 ,\ . . dx n , be reduced to the form 



(I.) 



*Wp» — Art C ^»+l~rAi2^n+2 • • • ""rAir^n+rj ^ 

the coefficients A n , A 12 , &c. being functions of the variables. 
Let P=c be an integral of the system. Then 
dP, , dV , , dP 



' CbX„±. v —\)» 



Substituting in this equation the values of dx x , dx 2 , . . . dx n given by (I.), we have 

/ dP , . dP , . dP , A dP\ , 



+ (—■ 

' \dx n+2 



dP 

*dx l 



■•^•i2^r+A, 



dP 
™dx a 



+A, 



dP 

n2 dx, 






+ 



/ dP 



A^4-A 



dP 

' dx 2 ' 



+ A -^;) 



dx n 



:0. 



As the differentials <fcr n+1 , <te re+2 , . . . dx n+r are now independent, we have, on equating 
their coefficients separately to 0, 



dP ^,A ^ p 



dxc," 






^p 



U/Xji 



+A -+A - +A ^-0 



i dx l 



% dx n 



(1X11 






dxc," 



+ A ^„-°> 



(II.) 



a system of r linear partial differential equations, the common integrals of which will be 
the integrals of the system (I.). We say ' the common integrals of which,' because in 
fact these equations express the conditions to all of which P must be subject in order 
that P—c may be an integral of the system (I.). 

The formal connexion between the systems (I.) and (II.) deserves to be carefully 
noticed. The several partial differential equations of the system (II.) may be formed by 
inspection from the columns in the right-hand member of the system (I.) by the follow- 
ing rule. For the differential dx n+i in any column write the differential coefficient 

dP . dP dP dP 

-^ — -, to this add the series of differential coefficients -r- > ^r, ... jr- multiplied in succes- 
sion by the descending coefficients of the column, and equate the final result to 0. 
The symmetrical form of the equation 

dJ, * r, +^* f « , --+si: <to »+' ==0 
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shows that by an exactly similar rule any system of n partial differential equations of 
the first order, the terms of which consist of the differential coefficients of P multiplied 
by functions of the independent variables # l5 # 2 , . . . #„+,, may be converted into an equi- 
valent system of r common differential equations of the first order. 

For the proposed system of partial differential equations is by algebraic reduction 
expressible in the form 



dP 



——A 



dP_ 
dx% 

d 2_ = a _ _ 



:A, 



11 da- , 

dP 

d%n+l 

dP 



4-A» 



+A 22 

■"-Hi 



dP 

d%n+Z 
dP 

dx n+2 ' 
dP_ 

+2 



dx n 



• +A lf 

• +A 2r 



dP 

dx„ +r '> 

dP 
dP 



(III.) 



If the values of ^ -^, ••- -^ in this system be substituted in the previous general 
equation, and the coefficients of the differential coefficients j — —, -i — 5 ••■■? — in the 

^ * ax n+1 «te n+2 dx n+r 

result be separately equated to 0, we shall have 

^'»+i"T" J ^-ii w^Pj-piiji 0M? 2 . . ,-y-A^i CW„— U, 

^^n+2 i ■"■12 OWiT-^ss ^*^2 • • •~\~-"-n2 ^^» == '-'5 r • » • * ■» • • • ^"*-*v 

^n+rT^lc 0/X\ -\-A-2r ("<£% • • •\-"-nr Ow M =U. 

These, equations may in like manner be formed by inspection from the columns of 

dP 
the second member of (III.), by writing for -* — in any column dx n+i , adding to this 

dxu dx 2 , . . . dx n multiplied by the descending coefficients of the column, and equating the 
final sum to 0. The rule for the one case differs from that for the other only in that 
differentials take the place of differential coefficients. 
As an objection may be felt to the legitimacy of that step of the above process in 

dP dP dP 
which, the differential coefficients ^-» ^- 5 ""fa" ^ em g eliminated, the coefficients of the 

remaining ones are separately equated to 0, I will point out another mode of procedure 
which leads to the same result, and which is founded upon Lagrange's method of 
solution. Let the equations of the system (III.) be added together after having been 
multiplied respectively by X„ X 2 , . . . K n , which are to be regarded as indeterminate func- 
tions of the variables w u # 2 , . . . x n+r . The result will be a linear partial differential equa- 
tion of the first order, of which the Lagrangian auxiliary system of ordinary differential 
equations will be 



d3C\ CUjOq 



&3Cn 



-dx n +i 



— dx, 



n+r 



A, Ajj Xj . . . + A„i A„ Ai r Aj . . . + A nr A„ 

Hence, eliminating \„ X 2 , . . . K n , or more strictly speaking the ratios 

Aj A 2 A n _i 

iT' I"' ' " T~' 

3o2 
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we have 

A n uX l .. . -j-.A. nl <&£„= — w3? m +d 



which agrees with the system (IV.). 

Of the Determination of the Number of Integrals of a system of Differential 
Equations of the First Order and Degree. 

We still suppose the given system of differential equations to be expressed under the 
general form (I.), and reduced by Prop. I. to the equivalent partial differential system (II.). 

Now if for the expression of that system we introduce a series of symbols A„ A 2 , . .. A, 
denned as follows, viz. 

A ^ I « ^ ■ A ^ . « ^ ,-, \ 

A '=S^+ A "^+ A --& 2 --+ A »^' ...... • (I-) 

the system will assume the form 

A 1 P=0, A 2 P=0,...A r P=0, (2.) 

and we shall now establish the following proposition 

Proposition H.—If AjPrsO, AjPrnO represent any two linear partial differential 
equations of the system (II.), then will the equation of which the symbolical expression is 

(AA-A,A)P=0 (3.) 

also be a linear partial differential equation of the first order, and it will be satisfied by 
all the common integrals of the two equations from which it is formed. 

For, representing any one of the quantities # 15 x 2 , . .. x n by x, and any function of those 

quantities by X, A t consists of a series of terms of the form X ~r. Again, representing 
any one of the same series of quantities #„ x 2 , . ...» n by y, and any function of them. by Y, 

A 

dy' 



Aj will consist of terms of the form Y -r-. Hence ( Aj A 2 — A 2 A,)P will consist of terms 



of the form 



[^dx * dy~ Y dy^dx)*' 



Effecting the differentiations, this term becomes 

X— — lYV— _ Y — ^ P — YX d * P 

dx dy ' dxdy dy dx dxdif 

or 

X— ^?_v— — 
dx dy dy dx 

which involves only the first differential coefficients of P. Hence (3.) will be a linear 
partial differential equation of the first order. 

Hence also the ultimate form of (3.) will be the same as if A j5 when operating on 
A,P, operated only on the coefficients A y . . . A nj involved in A,-, and vice versti. Thus the 
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ultimate form of (3.) will be 

dP dP dP 

(A ( Ay-A^«) £ j.+(A l A^-A,A - ) 3 ^ +(A ; A„,-- 4^)^=0. . , . (4.) 

Secondly, the equation (AAj— A J A,)P=0 will be satisfied by any common integral of 
A;P=:0 and A,P=0. 
For let <p =c be a common integral of the latter equations. Then, identically, 

A#=0, Ajp-0; 

therefore, since A,- and A ; involve only operations of differentiation together with alge- 
braic ones, 

A,A0=O, A,A^=0; 

whence (3.) is also identically satisfied. 

Proposition III. — If, by the successive application of Prop. II., and by permitted pro- 
cesses of algebraic elimination, we derive from the system of partial differential equations 

A,P=0, A 2 P=0, A,P=0, 

into which the given system of differential equations has been converted, a final system of 
partial differential equations which, while including the above system, shall be such that 
the application of Prop. II. to any pair of the equations contained shall lead only to an 
identity, then the number of integrals of the given system of differential equations will 
be equal to the number of variables they contain, diminished by the number of partial 
differential equations of the above final system. 
The developed form of the system 

A 1 P=0, A 2 P=0, ... A r P=0 , . (1.) 



is the following :- 



dP . . dP , . </P_ ft ^ 



dP . x dP . ^P_ ft 



#cfc* ?i .{-2 (tSO-t Oj95fi 



dP . . dP . . dP n 



(2.) 



Cw n ^.yf (IX 'j UltX/ft 

Comparing these with (4.), Prop. II., which is the developed form of the equation 
(AjA— A J A,)P=0, we see that the latter equation is necessarily algebraically independent 
of the above system ; for no equation derived from that system by algebraic processes 

dP dP 

could be free, as (4.) is, from all the differential coefficients -^ — ... -^ — . 

Again, as (A,A/— A 7 AJP=0 is satisfied by all the common integrals of A;P=0 and 
AyP=0, it follows that the system of r+1 equations, 

AJ^O, ... A,P=0, (A^,— AA) P =°> • (3-) 

will be satisfied by all the common integrals of the system (1.). To this system of r+1 
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equations we can also give a form analogous to the developed form of the system (2.). 

It will be noticed that the differential coefficients -= — .*. -= appear there, each in 

da:„+i dx n+r 

only one equation, and each with the coefficient unity. Now let the last equation of 
(3.) in its developed form be divided by the coefficient of -j-, and let also the value of 

dP 

•t— which it gives be substituted in the other equations of (3,) ; then we shall have in 

the whole a system of n-\-l equations possessing the same general character as the 
system (2.), To this new system the same procedure may be applied, viz. the genesis 
of a new equation by means of Prop. II., and the transference of another differential 

dP 

coefficient -r- — to the list of those which form the respective first terms of the equations 

of the system. We will suppose this procedure to have been repeated until a system 
composed of m partial differential equations such that the further application of Prop. II. 
leads only to identities has been formed. If n-\-r— m=p, that system will be of the 
form 



dxp+i ' u dx x ' pl dx p ' 

dP ■i-n,/ p ...+ii/ p -o 



dx p + m ~T~n- lmdx ^ ... -\-i± pmdw —v. 



(4.) 



And if, in analogy with former notation, we write 



d . TT d 



d^ i +B - u dJ+ 1I ' i d^ ••• + H ^ =A >-' 

it will take the symbolical form 

A,P=0, A a P=0, ...A m P=0; 

but it will differ from all former systems of equations in that all the conditions repre- 
sented by 

(AA-AA)P=0 (5.) 

Will be identically satisfied, — satisfied, in consequence not of any ascertained peculiarity 
of the integral P, but of the constitution of the system of symbols A 15 A 2 , . . .A m . 

The course of argument has shown that the common integrals of the system (4.) will 
be identical with those of the parent system (2.). Now we shall show that the exist- 
ence of the condition (5.) renders the integration theoretically possible ; that the system 
ofp ordinary differential equations into which, by Prop. II., the final system of partial 
differential equations (4.) is resolvable admits of exactly p integrals. As p=n-\-r— m, 
this is to say -that the number of integrals is equal to the number of original variables 
diminished by the number of final partial differential equations. 
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The proof of this will consist of two parts:— 1st. It will be shown that, if a system of 
p integrals exists, the conditions represented by (8.) will be identically satisfied. 

2ndly. It will be shown that, when the conditions represented by (5.) are identically 
satisfied, the solution either of the final system of partial differential equations (4.), or 
of the corresponding system of ordinary differential equations, by a system of p integrals 
is theoretically possible. 

It will follow from these conjoined, that the number of actually existing integrals is 
exactly^?. 

1st. The system of ordinary differential equations corresponding to (4.) may be ex- 
pressed in the form 

<te 1 =H 1I <Za? J1+1 +H 1 ,&r J , +1 ... +H lBl <k p+m , ) 

• • • • - . (6.) 

uX p —li pl a£ p+l -jrtLp i dX p+i ... -\-tx p ^aX p+m . J 

Now suppose this system to have p integrals. Then, by means of these, x lt x 2 , ... x p can 
be eliminated from the coefficients H n , &c. in the second members, which will thus 
become exact differentials of those functions of the variables x p+l ... x p+m which express 
the values of x u ... x p . Hence we shall have the system of conditions 

1^:^-1^^ ( 7 -) 

Jc representing any integer from 1 to p, and i, j any integers from 1 to m, and the 
bracketed symbols of differentiation referring to B. kj , TL ki as transformed. Hence, the 
unbracketed symbols referring to the prior state of the functions, we have 



( d \ d <fe, d dx p d 

dx p+i J dx p+ i'' dx p+{ dx\ ",' ~'dx p+i dx p 

1 l TT jL ITT J_ 

~dx p+i ' T - " dx 1 '" I"* 1 !* dx p 



=A, 



\dx p+J )- A i'> 



In the same way 

so that (7.) becomes 

AjH/y — A ? H A j=0. (8.) 

Now if we construct, in analogy with (4.), Prop. II., the developed form of the con- 
ditional equation (5.) of the present section, we shall have 

(A,H,- A,H„)~ . . . +(&M PJ - W g=0, 

or, 2 denoting summation from #=1 to k—p, 

2(A i H ¥7 A / H M )^ = 0. 

Now as by (8.) the coefficients vanish identically, the equation, and generally the system 
of conditional equations of which it is the type, will be identically satisfied. 
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2ndly. We proceed to show that the system of m linear partial differential equations 
(7.), represented under the form 

A l P=0,A 2 P=0, ...AJ?=0, 

and satisfying identically the system of conditions represented by 

(aa~aa)p=o, 

will admit of p integrals expressing distinct values of P ; and the system of ordinary dif- 
ferential equations (6.) corresponding to the above system of partial differential equations 
will be expressible as a system of exact differential equations, and will by integration 
give the above systems of integrals. 

Beginning with the first partial differential equation, of the system (4.), and forming 
the corresponding Lagrangian system of ordinary differential equations 

Cte p +1= H^-" = IL7 1 ' 

we see that the integrals of this system will be of the form 

w, ==<?! ...u p =c p , 

if, ...u p being functions of all the variables x v . . . x p+m , among which, by virtue of the 
integrals of the second line, x p+2 , . . . % p + m may be regarded as constant. The general 
integral will be 

J:(Mj, . . . Up, &p +i , . . . 0P p+m )—\), 

the form of F being perfectly arbitrary. 

Now the general form of any equation of the system (4.) is 

4^+H„f i ...+H„||=0 (9.) 

Let us transform this by assuming 

Wj . . . Up, CCp +1 . . . X p+m 

as independent variables. Then referring the right-hand members of the following 
equations to the new, the left-hand to the old system of variables, we have 

d? _ <?P d? <foi_ d? du p 

dxp+i dx p+i ~1~ du x dxp+i " ' ' du p dx p+ i 

— — ^l^h i §1 d Jh 

d%i ~ du x dx x ' ' ' "r du p dxj 



dxp du x dx p •'•"•" du p dx p ' 
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whence, substituting in (9.), 



dx p +i * du 



7-A l TT *fl .XT ^l\ 



_1_— /_^__I_TT ^f? J_TT ^A ft, 

"T"^ \& p+i -r **■* dx l • ' ' + n * «fcj — u ' 
or 

JP I /A \ rfP , /A ^ P A /IAN 

^+( A ^)^--.+(A^)^=0; (10.) 

and in this equation we may give to i the successive values 1,2, ...p. If, then, we write 

4^ +(AiWi) i • • • +(***) x-=x* 

we see that the proposed transformation will have the effect of converting 

A„ A 2 , ... A p 
into 

a;, a;, ... a; 

respectively, and the system of partial differential equations into 

a;p=o, a;p=o, ... a;p=o. 

The developed form of the first of these equations is 

But since w„ ... u p are integrals of AiP=0, we have 

A,w,=0, ... A,^— 0, 
so that the equation AiP=0 reduces to 

dx pJrl 

"We learn from this that x p ^ will not explicitly appear in P after the transformation 
which introduces w x , ... u p . 

The developed form of the remaining^— 1 equations represented by (10.) will be 



dP i fK ^ p , /a ^ p A 
^+(A 3 <^-...+(A 2 ^=0, 



-p 






(ii.) 



and we shall next show that the variable x p+1 will not present itself in the coefficients 
(A 2 «,), &c. 

The general form of such coefficients is A t u s , where i has any value from 2 to m, and 
s any value from 1 to p. 

Now if in (5.), which is true independently of the nature of the function P, we make 
mdccclxii. 3 p 
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j=l and P=w s , we have 

(AA-AAK=0. 

But AiW,=0; therefore, by the above, 

A,AA=0. 

Now AjW s is expressible at most as a function of m„ ... u p , s p+1 ... x p+m . But this trans- 
formation converts, as has been seen, A, into -j- — . Thus we have 

-A s w s =0, 



dxp +1 

so that A t u s is free from x p+1 . Thus the system (11.) is free from # p+1 . 

Lastly, since by the above transformation A 2 ... A p are converted into A 2 ... A^, the 
system of conditions (AA - AyAj)P=0 is converted into (A^A]— AA)P=0. 

It is thus seen that the system ofp partial differential equations 

A,P=0, A 2 P=0, ... A p P=0, 
containing p-\-m independent variables a?„ # 2 , ...%. m , and satisfying the conditions 

(AA,-AA)P=0, 
is convertible into a system of p— 1 partial differential equations, 

a;p=o, a;p=o, a;p=o, 

containing p+m— 1 independent variables w, ... w p , # p + 2 ... a^ +m , and satisfying the con- 
dition 

(aa-aa)p=o. 

And as this system possesses the same character as that upon which the previous trans- 
formation depended, it will admit of transformation into a system ofp— 2 partial differ- 
ential equations containing p-\-m— 2 independent variables; and so on until we arrive 
at a single final partial differential equation containing p+1 independent variables, and 
having therefore p distinct integrals, which will be the common integrals of the primary 
system of partial differential equations as well as of the system of ordinary differential 
equations to which they correspond. •"=.--.'-' 

Cor. The property of the coefficients A 2 «„ &c. of the system (11.), of being free 
from the variable w p+1 , enables us, by properly determining the integrals of the par- 
tial differential equation A,P=0, to reduce the system to a form of great simplicity. 

Let A t Uj be any one* of those coefficients. • Its developed form is 

(^+ H »d^"-+ H *Er)% *•••■-.. (12.) 

Now as this expression will, after the performance of the differentiations, be free from 
#p+i> and as the differentiations are none of them with respect to ce p+1 , we can give to 
x p+1 in it any particular value before differentiation without affecting the final result. 
Let us then suppose that in H H , ... H^, and in u p w p+1 is made equal to 0. Now it is 
possible so to determine the integrals u f as functions of the variables x ls ... # p+r , that 
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when # p+ j==0 each % shall reduce to x f . For this purpose it is only necessary to choose 

as arbitrary constants 'a set of arbitrary values of $„ x 2 , ... x p corresponding to x p+1 =0. 

]Let w"i, x' 2 , ... ~x'p be such arbitrary constants, and let the given system of integrals be 

reduced to the form 

u % =x[, ... u p =x' p , 

and the functions w n ...'% will possess the required property. Changing, then, each 
Uj into Xj, the expression (1,2.) reduces to H ;i , and it only remains to express this 
in terms of u x , ...%; which, as ^ +1 =0, is done by merely changing Xi,„X p into 

U t , :.-.:. U pt 

Thus the system (11.) is reduced to 

dP . /TT ' dP . m , rfP 



&^ +(H) ^ s; ••• +(Hp2) ^ =0, 



+( H »»te ••• +( H ^)& = ' 



Where the brackets denote that in the enclosed portion x p+1 is to be made 0, and 
Xx, x 2 , ... ^converted into u„ u 2 , .... w p . 

Now this form is identical, the above conversion of letters excepted, with that of the 
system (4.), omitting that equation of the latter system by the integration of which the 
forms of u x , r ...u p are determined. 

It follows from the above that, obtaining the integrals of A,P=0 in such a form that 
the arbitrary constants shall represent the arbitrary values of x x ■ .. x p when x p+l =0, and 
representing the functions which are equal to those arbitrary constants by x[,x a , ...x' p , 

then if in the remaining equations t A 2 P = ... AJP=0 we change .x l ,.. : ..x pi _-,«..— to 

(tSG u90p 

■ d ' d ' dP 

x\y...x p , jt, • • • -r-r, and x p+l and' r- — to 0, the common integrals of the transformed 

dX\ dx p "*p+i 

system of p— 1 will be the f same as those of the previous system of p partial differential 
equations. In the same way a third system of^> — 2 partial differential equations may 
be formed, and so on, till we obtain a single final partial differential equation which will 
have the common integrals of the parent system. By this method, which is due to 
Jacobi and Nutani, all the labour of the successive transformations is avoided. The 
successive integrals thus introduced are termed ' Haupt-integrale.' 

Instead of applying the foregoing methods,, general or particular, to the final system 
of partial differential equations, we may apply it to the solution of the corresponding 
final system of ordinary differential equations. In this case they would really represent 
the method of solution known as the variation of parameters, and the conditions 
(A^A,— A ; A^)P=0 would secure the sufficiency of that method. If in the system of 
ordinary differential equations, (6.) we regard # p+a , .. . x p + m . as constant, we get 

dXi—B.u^Xp^ — 0^ 

(13.) 

dx p — H n dx pJr i = 0. J 
3p2 
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Integrate this in the form 

«,=<?„ ...u p =c P ; 

then, treating (?„... c p as functions of the variables before regarded as constant, and 
endeavouring to satisfy the unreduced equations (6.), we obtain, in virtue of the above 
conditions, a system of differential equations equal in number to the system given, but 
containing one variable fewer. The system (13.) by which the forms of m„...m p are 
here determined, is the Lagrangian auxiliary of the first partial differential equation 
A 1 V=0 integrated in the other method; and so in each subsequent stage. And with 
respect to the other parts of the process, it obviously makes no difference whether we 
take as the new variables «„ ...u or c lt . ... c p under the condition (necessarily involved in 
the method of the variation of parameters) that they shall after integration be replaced 
by «„ . . . u p . But it would not have sufficed simply to refer the solution of the final 
system of ordinary differential equations to the method of the variation of parameters, 
first, because the necessity and sufficiency of the conditions which form the ground of 
that method and are the warrant of its success were to be shown ; secondly, because the 
Connexion of the systems of ordinary differential equations which arise in the method of 
parameters with the successive partial differential equations forms an essential part of 
the demonstration. 

General Bute. 

The results of the foregoing inquiry may be collected into the following General 
Rule : — 

To find the number^? of possible integrals of a system of n differential equations of the 
first order connecting n-\-r variables a?„ x 2 , ... x n+r , and to determine those integrals. 

ifo«fe.- i -Suppose P an integral of the given system. Determine from the given system 
dx r , dx s ,...dx n as functions of the other differentials. Substitute these values in the 
equation 

and equate to the coefficients of the remaining differentials. This will give a system 
of r partial differential equations of the form 

dP . . dP . dP 

+ A n^7. •••+ A i»^r— u > 



dx n+1 ~^ n dx x ~^ ln dx n 






(XiL j^-j-)- (XX \ (litXjyi 

r of the differential coefficients appearing each only in one equation and with coefficient 
unity. Representing these equations in the symbolical form 

A^O, A 2 P=0,...A r P=0, 

deduce any equation or equations of the form 

(AAj- 4AJP=0. 
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If all such prove to be identities, the given system of differential equations admits of n 
integrals, and is reducible to a system of exact differential equations. But if any such 
equation is not an identity, it will constitute a new partial differential equation of the 
form 

B — 4-B — 4-"R —-0 

And this, combined with the previous ones, will enable us to form a system of r+1 

dP dP dV 

partial differential equations, in which -5— 5 -7 > ... -. — ■ appear each in only one equa- 

tion and with coefficient unity. Upon this system let the same process be repeated as 
upon the previous system of r partial differential equations, and so continually repeated 
until we arrive at a final system of partial differential equations such that, if that system 
be represented in the form 

A 1 P=0...A m P=0, 
the condition 

(aa-aa)p=o 

shall be identically satisfied for every pair. 

Then, the number of such partial differential equations being m, the number of inte- 
grals of the original system of partial differential equations will be n-\-r— m, i. e. it will 
be equal to the number of the original variables diminished by the number of final 
partial differential equations. 

And if by that final system we eliminate m of the differential coefficients from 

and equate to the coefficients of the remaining differential coefficients, we shall have 
a system of %-\-r— m differential equations expressible as exact differential equations for 
the determination of the integrals. 

Actually to determine these, we should endeavour in the first instance to reduce the 
final system of differential equations, as such reduction is theoretically possible, to a system 
of exact differential equations. If the means of doing this are not obvious, the method 
of the variation of parameters or the equivalent methods of Prop. III. must be applied. 

Lastly, if the process which consists in the application of the theorem (A^A,— A j A i )P=0 
do not stop with the formation of the final system of partial differential equations, but 
lead to algebraic relations among the variables, the given system of differential equa- 
tions will have no integrals properly so called, but it may admit of solutions analogous 
to those the theory of which has been developed by Ppaff, Jacobi, and others for the 

differential equation 

Xj^i+X^a^ . . . +X a $r a =0. 
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Applications. 

1st. Suppose it required to find the number of integrals of the form P=0, which the 
system of differential equations 

dz—(t +wy-\-wz)dw-}-(xzt-\-y—wy)dy, 

dt=(y+z — Bx)dx-\-(zt —y)dy 

admits, and to determine such integrals. 
Eliminating dz and dt from the equation 

dP , ,iP.' dP , dP . n 
-djdw+^dy+^dz+^dt^O, 

and equating to the coefficients of dz and dt in the result, we have 

— + (t+ X y+xz)-^+(y+z-Zx) Tt =0, ... . . . (1.) 



Hence writing 



f+M+^-^)S+(^-^)f=° ( 2 



&'+~+(wzt+y-wy) j z +(zt-y) J, 

and forming the equation 

(AA'-A'A)P=G, 

we have, on rejecting a common algebraic factor, 

dP . dP n ... 

^^+^=°- • • ' • ( 3 

dP 
By substituting in (1;) and (2.) the value of -^ hence obtained, we have the system of 

three equations, 

dP dP =Q 
dy " dz ' 

dV dP 

-dt +W dz- = °- 

jNbw if upon any two of these we repeat the same process as upon (1.) and (2.), we 

obtain as the result 0==0. Thus the system of partial differential equations is complete. 

As then there are three equations in this final system, while the number of original 

variables was four, the primary system will admit of one integral of the form P=c, 

dP dP dP 
To obtain this integral, eliminate -^p -±-, -^ from the above equations and 

^dx+ d ^dy+~dz+~dt=0, 
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and equate to the coefficient of -7- in the result. We find 

dz—(t+ 3a?)d$—gdy—a;dtY=0, 
the integral of which is 

z— xt— x 3 — y*=c. 

2nd. The solution of the partial differential equation 

BH-Ss+Ttf+U(s 2 -rt)=V, 

as well as of the special equations 

Er+Ss+Ttf=V, 

the theory of which constitutes an exception to that of the more general form, depends, 
in general upon the integration of three simultaneous differential equations between 
five -variables. ' To this"integration the method of the foregoing sections is applicable. 

The only cases, for which the .theory of the ultimate solution can be said to be com- 
plete, are those in which the auxiliary system of common differential equations admits 
either three integrals of the form P==<?, or two integrals of that form. 

We may apply the method of the foregoing sections, not only to the determination 
of the integrals, but also to the discovery of the a priori conditions connecting the coefEU; 
cients K, S, T, II, Y in order that each of these species of integration may be possible. 

For example, the solution of the equation 

Er+Ss+T£+(s 2 -f£)=V 
depending upon the integration of the system 

dq=— m-fix+lkdy, , , „ 

dp= — m 2 dy-{-1dx, 

dz = pdz -\-qdy, 
in which m^ and m 2 are roots of the equation 

m 2 -Sm+KT-V=0, 

let it be required to determine the conditions under which the system admits three 
integrals. 

Eliminating dq, dp, dz between the above equations and 

dV, , J d? , .iS , . <ZP, , iP, . 

and equating to the coefficients of dx and dy in the result, we obtain two partial 
differential equations which may be thus represented, viz. 

AP=0, A'P=0, 
in which 

d d „, d d 

., d TO d d d 
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That there may be three integrals, it is here necessary that there should be but two 
partial differential equations in the completed system. Hence the equation 

(AA'-A'A)P=0 
must vanish identically. 

Developing this, we have the conditions 

AE+A'm,=0, Am 2 +A'T=0, A 2 -A'_p=0. 

Now on performing the operations denoted by A and A', the last equation gives 

m 3 — m,=0. 
Hence referring to the quadratic, we see that 

S 2 -4RT+4V=0 (I.) 

To this must be added the two other reduced conditions, 

AR+A'm=0, (II.) 

Am+A'T=0, (III.) 

m representing one of the equal roots of the reduced quadratic. 

The first of the above conditions was given by Ampere*. The others are probably 
new. Satisfied, they enable us to predict that the partial differential equation under 
consideration admits a complete primitive involving three constants, and a general 
primitive arising from the variation of those constants in subjection to any two arbitrary 
conditions. 

3rd. "We have supposed each linear partial differential equation employed in the pro- 
cesses of this paper to be of the form 

, rfP Li dP , . <*P _ ft 

Ai ^; +A *dx 2 - - - "t A » dx n - u ' 

and we have supposed each system of partial differential equations which arises, to be so 
reduced that each equation shall have some one of the partial differential coefficients of 
P entering into it alone and with a coefficient equal to unity. 

The first of these conditions is virtually sufficiently general, because any linear par- 
tial differential equation can be deprived of its second member. The advantage of the 
second condition is that each newly-formed equation will be really new, and not an 
algebraic combination of the old ones. 

But neither of these conditions is necessary. From two linear partial differential 
equations of the form 

A 1 P=H, A 2 P=K, 

in which H and K are functions of the independent variables, arises a new equation, 

(A 1 A 2 -A 2 A 1 )P=A 1 K-A 2 H, (1.) 

which will be satisfied by all the simultaneous integrals of the equations from which it 
is derived. 

It may be rigorously proved that, in applying this process, the generated system 

* Journal de l'Ecole Polyteclmique, Cahier 18 e . 
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(including the original equations) will be complete when no new equation arises from 
the combination of any one of the equations with any one of the equations ©f the 
original system. 

I will illustrate this by investigating the conditions of integrability of the expression 



Vy*' y> dx> dx» '" dx") aw ' 



If this expression admit of an integral V, it is easy to see that V will satisfy the two 
partial differential equations 

dV , dV , dV , dV 

&+^+^».+y»^==?>> (2-) 

£ =°' <*> 

in which 

dy dry 
y x ~Tx> y*~~d1tf- ' 

and <p stands for 

<p{x,y,y„...y n ). 

The above are, in fact, the partial differential equations which we should obtain by 
Prop. I. as the equivalents of the system of ordinary differential equations, 

dV=<pdx, 

dy=y t dx, dy l =y 2 dx, ... dy n ^=y n dx. 

If we write 

(d\ d d d d 

Tx) —Tx+yHy+y^ '" +#»"^I? 

the above partial differential equations become 

(s)v=*...(L), £=0...(II.) 

The combination of (I.) with (II.) (by the theorem (1.)), then of (I.) with the result, and 
so on, gives a series of equations which may be thus expressed : — 

£r=^ (^) 

£t=^_a . (iv.) 

f = A ^ (V - } 

= A <p, (VI.) 

in which 

A d__ /d\ d /dy d _ &c 

r ~~dy r \dx)dy r+l '\dx)dy r+2 *' 

dy \dxjdy 1 "^ydx} dy% 
MDCCCLXII. 3 Q 
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The combination of (IT.) with (III.) ... (V.) gives the series of conditions 

The conditions of integrability are expressed by (VI.) and (VII.). These satisfied, the 
equations (III.), (IV.), ... (V.) show that <pdw can be expressed as an exact differential 
with respect to cc, y, y„ y n ^, in the form 

$dx={$—yA$—yA& ... —ynK<p)&c 

+A 1 <pdy+A !i <pdy 1 ... + A n <pdy n _„ 
a result first established by M. Sakeus. 



